Interaction of spatial solitons in nonlinear optical medium 
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The effects caused by nonresonant nonfinear interaction between noncoilinear seif-focusing beams 
are considered in 2D opticai sampies using multi-scale analysis. The analytical expression for the 
beams trajectories shift due to the mutual interaction is derived, and the range of parameters is given 
beyond which the mentioned consideration fails. We compare our results with the naive geometrical 
optics model. It is shown that these two approaches give the same results. This justifies the use of 
geometrical optics approach for description of elastic and almost-elastic collision processes both in 
Kerr and saturable nonlinear media. 
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Nonlinear localized waves, or soliton-like excitations 
play important role in many branches of physics: non- 
linear optics 0|, plasma physics 0, hydrodynamics 0, 
magnetic systems etc. In contrast with linear excita- 
tions such nonlinear creations may be exceedingly stable. 
That is, they can propagate over long distances with- 
out distortion. However most exciting feature of soli- 
ton phenomena is their interaction processes. In partic- 
ular, when they collide with other ones, solitons exhibit 
partical-like behavior ||. 

Despite the fact that there is a large diversity of non- 
linear physical systems exhibiting soliton like excitations, 
due to universal properties of such creations, nonlinear 
localization dynamics can be described only within a few 
theoretical models |2(. This fact is of great importance. 
In particular, this allows one to study experimentally 
nonlinear phenomena in most convenient physical sys- 
tems, while the direct experimental investigation of the 
particular system the subject of interest may be more dif- 
ficult or even impossible. To date, such "model" experi- 
mental systems are nonlinear spin waves in ferromagnetic 
films 3 and spatial optical solitons ||. However, mag- 
netic envelope solitons can be observed only in (quasi) 
one dimensional samples and due to transverse instabil- 
ities they are unstable in higher space dimensions 0. 
Thus recently suggested interaction effects of non- 
collinearly propagating ID envelope solitons in 2D mag- 
netically ordered samples doubtfully will be easily real- 
ized experimentally. On the other hand nonlinear opti- 
cal medium is most appropriate for the mentioned pur- 
pose (study ID solitons noncoilinear interaction in 2D 
samples). In particular, spatial optical solitons exhibit 
richness of characteristics not found for localized waves 
in other nonlinear media [jj. Indeed, much theoretical 
and experimental investigations have been performed for 
optical spatial solitons: elastic interaction of Kerr soli- 
tons 0, H, IS EH , almost elastic and inelastic collisions 
of solitons in saturable media including fusion, fission, 



annihilation, and spiraling occurances 0] (see also 0). 

In the present paper we consider theoretically the prob- 
lem of nonresonant interaction of Kerr spatial optical soli- 
tons. The method used here has been suggested for ana- 
lytical description of interaction of noncollinearly propa- 
gating ID envelope solitons of magnetization in 2D mag- 
netically ordered samples . This method itself is a gen- 
eralization of well-known ID multiple scale perturbation 
theory 0,0] for higher space dimensions. Since this ap- 
proach allows us to study the case of interaction of two 
spatial solitons with different carrier wave frequencies, 
for the nonresonant interaction of solitons the results pre- 
sented here are more general compared to ones obtained 
in 0, where the exact solutions are found but they con- 
cern only the case of spatial solitons interaction with the 
same carrier wave number. Later these analytical re- 
sults has been used for suggesting different applications 
of spatial soliton interactions, e.g. producing nonlinear 
photonic switching and all optical logic elements . 
Optical solitons dragging and collisions in the presence of 
absorptions has been also studied 0] applying numeri- 
cal methods. In the present paper we obtain analytical 
results for collisions of spatial solitons with different car- 
rier wave number and point out possible relevance of this 
study to the above-mentioned applications. 

For simplicity we consider the interaction of spatial 
optical solitons in isotropic thin optical films and assume 
that the electric field is normal to the film plane. If the 
medium is off-resonant with respect to the optical sig- 
nal and the optical film is thin enough, then the dis- 
persion can be neglected and modulations develop only 
along (single) transverse direction 0,0|. As a result, 
in the case of appropriate sign of the nonlinear coeffi- 
cient, so called ID spatial solitons (self-focusing beams) 
are formed in 2D samples. Obviously one can consider 
the crossing of two beams and study analytically the in- 
fluence of one self-focusing beam on the other one using 
the above mentioned method 61. 
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One could try to understand the nonlinear effect of spa- 
tial solitons interaction from the naive geometrical optics 
model: the propagation of the intense beam through the 
sample locally causes the increasing of the refraction in- 
dex due to the nonlinear reaction of the medium (Kerr 
effect). This in turn generates wave-guiding area and as 
a consequence the beam becomes self- focusing [jj. At 
the same time, as long as the refraction index within 
the beam area is bigger than outside it is natural to ex- 
pect that the second beam will be bent during crossing 
the first beam area and eventually its trajectory will be 
shifted after the interaction as is shown schematically in 
Fig. 1. From the same geometrical optics consideration 
it follows that this shift should be zero for perpendicular 
beams. However, the interaction process is much more 
complicated. Actually the second beam affects the in- 
duced waveguide of the first beam (the first (wide) beam 
in Fig. 1 is slightly shifted as well). This gives rise to 
" self-action" of the beam through other one during the in- 
teraction. In addition, interference effects may take place 
between the carrier waves of the interacting beams. For 
Kerr spatial solitons it is well established that for large 
enough converging input angles the solitons pass through 
each other unaffected. The only effect of such nonreso- 
nant interaction is the trajectory shift of the interacting 
beams j^J- Thus, qualitatively, the effect is the same as 
it comes out from the naive geometrical optics consider- 
ation. Surprisingly, in this paper we find that the results 
are the same even quantitatively. 



linearly on electric field E as follows: 



(1) 



where x*- 1 ^ and x^ are linear and nonlinear polarisabil- 
ity constants respectively. For simplicity we consider 
here centrally symmetric materials and therefore from 
the symmetry reasons the second order term in Exp. 
is identically zero. The wave equation for the electric 
field reads (see for more details e.g. 0]): 



W 2 E - pE tt = 7 (S 3 ) it 



(2) 



with coefficients /3 = (1 + Attx^) / c 2 and 7 = Airx^ /c 2 . 
The nabla operator acts in 2D space as long as film sam- 
ples are considered in the present paper. 

Let us consider the weakly nonlinear case i.e. when the 
cubic term is much smaller than the linear one. We do 
not repeat all the steps of calculations to obtain the ID 
spatial soliton solution, let us just mention that a weakly 
nonlinear wave with a slowly varying envelope is sought 
in the following way [12| : 



a=l 



il{kr—ujt) 



(3) 



l= — oo 



where frequency to and wave vector k of carrier wave are 
connected via simple dispersion relation u> = k/^/]3; the 

envelope <pj^ — (tp^] is a function of slow variables 



In nonlinear Kerr medium polarization P depends non- £ = e{f — vt) and f = e 2 f/2k; v = du/dk is a group 

velocity of linear wave v \\ k and e is a formal parameter 
defining the smallness or " slowness" of the physical quan- 
tity before which it appears. Then building the ordinary 
perturbation scenario [2j, [12( in the third approximation 
over e one gets ID nonlinear Schrodinger (NLS) equation: 







n 2 







dr dS, 2 



(i) 
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= 0, 



(4) 



FIG. 1: Schematic picture of the interaction process between 
self-focusing beams in off-resonant optical medium. Solid 
lines indicate "borders" of the beams, qo defines an angle 
between first (narrow) beam and normal vector 712 of the sec- 
ond (wide) beam; k\ and ki are their carrier wave vectors, 
respectively; 5li stands for a shift of trajectory of the first 
beam, which is caused by the nonlinear interaction effects. 
Let us mention that the second beam trajectory is also slightly 
shifted. 



which has well-known spatial soliton (self-focusing beam) 
solution if 7 > 0. Physically, the spatial soliton forma- 
tion is the result of balance between defocusing diffrac- 
tion and focusing nonlinearity. It is worth to note that, 
since diffraction in general is strong, the nonlinearities 
involved in spatial solitons are much stronger compared 
to nonlinearities involved in temporal solitons. As we 
see from @ the soliton envelope in the leading approx- 
imation is a function of variables £ = £j_ = e(nr) and 
t = fji = e 2 (kr)/2k 2 where r plays a role of "time" and 
ft is a unit vector perpendicular to k. Let us empha- 
size once again that we have got ID NLS due to fact 
that the longitudinal dispersion in off-resonant optical 
medium could be neglected, in other words wave group 
velocity v = duj/dk does not depend on wave number k. 
In case of presence of longitudinal dispersion the physical 
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process would be described 2] by 2D NLS for which ID 
solitonic solution would be unstable. 

Now we shall start a main task of the paper, particu- 
larly, the analytical investigation of interaction between 
noncollinear self-focusing beams. For that purpose we 
are seeking for the solution following the general method 
developed in Refs. 0: 



E 



E 

hh=—i 



e^'i^-^i^+^iO^ (5) 



where uii ± i 2 = hui + 1 2 uj 2 ; fa 1 i 2 = hk\ + fak 2 ; £*>i, fa 
and W2, k 2 are carrier frequencies and wave vectors, re- 
spectively; (ff"i and Q^iz are functions of slow variables 
(P=l,2) 



Hp 



(n p r) - eipp(£,i,£, 2 ,Ti,T 2 ) 



2fc2 ; 



(6) 



where n\ and n 2 are unit vectors perpendicular to carrier 
wave vectors k\ and k 2 , respectively. Proceeding further 
with the similar calculations as it was done in Refs. 
we come in the leading approximation to the following 
solution: 



E 



(i) i 



fel f—UJl 



t +n«) 



(1) i 
tin e 



k 2 f— u^i+ni 



+ c.c. 



{ (7) 

where "c.c" denotes complex conjugated terms; (p^f and 
tp$ are the solutions of ID NLS (see Eq. (gj) 



with 

derivatives over set of slow variables £i, n and ^2, i~ 2 , 
respectively. For example, in the simplest case of one- 
soliton envelopes ip^J and we have: 



<p§\ = \A 1 \-aech{\A 1 \ > /6^[(jt 1 r f )-4 1) ]}, 

{lAal^J^rO-^ 15 ]}, (8) 



\Ai\ ■ sechi 



and A\ and A 2 are amplitudes of the first and the second 
self-focusing beams, respectively. Besides that, phase 
and position shifts of the first self focusing beam induced 
by weakly nonlinear interaction with the second beam 
are defined by the following formulas: 



_ (gigg) ggio 
#6 (kin 2 ) 9 & 



3 7 ^|^| 2 



(nin 2 ) 
(kin 2 ) 2 



(9) 



As far as according to perturbative approach we have a 
following scaling dip^ /d£ 2 ~ e 1 and taking into consid- 
eration the dispersion relation (uj\ = k\j f3) the following 
restriction upon the soliton parameters is derived: 




(3 cos 2 ao 



«1, 



(10) 



where ao = (tt/2) — 6q, and 9q is an angle between the 
self-focusing beams. From Exp. © it is easy to get 
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FIG. 2: Optical ray refraction through the self-focusing beam 
area with borders denoted by horizontal solid lines. 



analytical expression for trajectory shift of the first beam 
caused by the nonresonant interaction with other one: 



Sh = 



4 1} (oo) 

37 sinao 
cos 2 ao 



■00 



(11) 



Particularly, in case of one-solitonic envelopes © one 
gets from (llljl the analytical expression for the trajectory 
shift of the first spatial soliton: 



Sh = \A 2 \ 



67 sinao 
f3u>2 cos 2 ao 



(12) 



In addition we want to emphasize that beyond the limit 
given by (|10fl the dynamics is governed by a set of two 
coupled NLS type equations, which in general is not in- 
tegrable model and gives qualitatively different behavior 
of interacting solitons (see 0] and discussion there). 

Now let us compare the results obtained above with 
the picture given by the naive geometric optics consider- 
ation. This will give us more deep insight to the problem. 
First suppose that one has only one self-focusing beam 
(particularly, the second (wide) beam) and let us calcu- 
late how it changes refraction index (see Fig. 2). In view 
of dependence of polarization upon the electric field 
we can write down the expression for refraction index as 
(j = yj3 + 3jE 2 , where E 2 denotes electric field in the 
self-focusing beam area and it has a form of envelope 
spatial soliton (jHJ. Averaging refraction index over fast 
variables r and t in the weakly nonlinear limit (the term 
proportional to E\ is small) we get the following approx- 
imate formula for slowly varying (along axis x) averaged 
refraction index u{x) = y/(3[l + {3j/2/3)eJ\ . 

Let us now solve the following geometrical optics prob- 
lem, particularly, how the optical rays refract propagat- 
ing through the area of the second beam with slowly 
changing refraction index. For that purpose we note that 
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an angle a between the ray and normal vector (with re- 
spect to the beam) at any point could be calculated via 
simple refraction formula: a(x)sma = uosinao, where 
do = vp- Taking into account the fact that the tra- 
jectory shift of the ray could be calculated as follows: 

oo 

Sli = cosao / dcc(tanao — tana), we come exactly to 

— oo 

the formula which we have obtained from multi-scale 
analysis 0. However, to avoid any misunderstanding it 
has to be stressed that the geometrical optics approach 
to the spatial solitons interaction problem is not self- 
consistent. Indeed, as we have pointed out already, in this 
model the " self-action" effects of the first soliton through 
other one are neglected. Physically this means that the 
first beam is linear. But the diffraction in the nonlin- 
ear problem considered here in not negligible. Thus the 
first beam will diffract and the interaction picture given 
by the geometrical optics model is not meaningful under 
realistic experimental situations of the spatial solitons 
interaction. 

Although the geometrical optics approach in some par- 
ticular cases gives full understanding of the solitons in- 
teraction processes fl8], in general one expects that this 
approach is valid only for qualitative description of the in- 
coherent spatial solitons interaction y|. As is mentioned 
above the nonlinear "self-action" of the beam through 
other one is neglected without justification in the geo- 
metrical optics model of solitons collisions. However, the 
analysis presented here shows that this additional nonlin- 
ear "self-action" during the interaction process does not 
affect the soliton dynamics asymptotically. This is why 
such naive geometrical optics model gives correct results 
even for almost-elastic collision processes between soli- 
tons in saturable media 

In summary, in the present paper we have considered 
the problem of nonresonant interaction of Kerr spatial 
solitons (sclf-focusing beams) with different carrier wave 
frequencies using multi-scale analysis. It is shown that 
the beams trajectories are shifted due to mutual interac- 
tion. The analytical expressions for these shifts are ob- 
tained as well. Surprisingly the naive geometrical optics 
model of the solitons collision is in full agreement with 
the results of general theory. In particular, this shows 
that the "self-action" of the soliton caused by nonreso- 
nant interaction process does not changes it's asymptot- 
ical behavior after the collision. This in turn justifies 
the use of the geometrical optics model for description of 
elastic and almost-elastic collision processes both in Kerr 
and saturable media. 
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